REMARKS 



Examiner 's Objections to the Drawings; Amendments to the Specification: 

Examiner has objected to the drawings under 37 CFR 1.83(a) as not showing every 
feature of the invention specified in the claims. In response thereto, Applicant respectfully 
directs Examiner's attention to the new drawing sheet submitted herewith. Applicant 
respectfully believes that the newly submitted drawing addresses Examiner's basis of the 
objection under 37 CFR 1.83(a) and, accordingly, requests Examiner's approval of same. 
However, should Applicant's drawing not be fully responsive to Examiner's objections and 
requirements, then Applicant respectfully requests a non-final opportunity to correct same. 
Applicant further submits that the specification has now been amended, as detailed hereinabove, 
to appropriately reference and described Applicant's new drawing. Applicant notes that no new 
matter has been presented in either the drawing or amendments. 

Examiner's Objections & Rejections to the Specification & Claims; Applicant's Responses: 

Examiner has objected to the specification under 35 USC 112, first paragraph; and, has 
rejected elected Claims 1-10, 13, 14, 89, and 91-95 under 35 USC 112, first paragraph; 35 USC 
101; 35 USC 112, second paragraph; 35 USC 103(a); and, 35 U.S.C. 102(b), citing a multitude 
of published articles and a series of patent references in support of Examiner's various arguments 
asserting that there is no reputable evidence of record to support any allegations or claims of 
invention made by Applicant, and that Applicant's invention is based on the concept of cold 
fusion. (Although Examiner has raised many additional related arguments, for purposes of 
brevity Applicant will not reiterate same, but will rather draw reference to general arguments 
reiterated by Examiner throughout Examiner's nearly fifty (50) page Office Action.) 
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In response to Examiner's arguments, and with all due respect to Examiner, Applicant 
respectfully submits that Examiner has largely misread, misrepresented, and misunderstood 
Applicant's invention, both in the specification and in the claimed subject matter. Indeed, 
although Applicant believes that Examiner has provided a lengthy and thoughtful response, 
Examiner's response is highly prejudicial and, for the most part, inaccurately represents the 
scope and nature of Applicants proffered technological and inventive advances. 

To begin, Examiner has given a lengthy history of cold fusion and the position of the 
scientific community in response thereto. Unfortunately, Examiner has improperly and 
inaccurately categorized Applicant's invention as that of a cold fusion process. Quite to 
contrary, however, Applicant's invention is a variant of hot fusion - a very different process 
from cold fusion, and recognized in the scientific community. For example, when electrons in a 
TV picture tube impact the screen exciting phosphors to produce light, one may touch the screen 
without burning one's hand, yet every competent physicist can calculate that the kinetic 
temperature of such electrons exceeds ten million degrees Centigrade. Applicant's invention is 
not a cold fusion process, nor was it ever claimed to be such a process. Applicant addressed the 
cold fusion process in his application merely to provide a historical synopsis of the various 
fusion reaction efforts (this is supported by the fact that Applicant presents the cold fusion 
description in the "Background of the Invention" section of the patent application - the same 
section from which Examiner draws support for his rejections). However, Applicant has not 
stated or claimed anywhere in his patent application that his invention is one of cold fusion, nor 
has Examiner been able to cite to any such claim made by Applicant. 
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Instead, based solely on a single paragraph contained on pages 7-8 of Applicant's 
"Background of the Invention" section (many more pages of which contain far more paragraphs 
describing fusion reactions unrelated to cold fusion), Examiner has elected to classify, much to 
the detriment and prejudice of Applicant, the instant invention as one of cold fusion, and has thus 
provided Applicant with an Office Action entirely based on such unfounded assertions. 
Additionally, the numerous references cited by Examiner, although applicable to cold fusion 
processes, are wholly unrelated and inapplicable to Applicant's fusion process. Applicant cannot 
accurately and properly respond to an action comprising seemingly endless arguments rejecting 
an invention that Applicant has not described or claimed as being one of his own. 

Scope and Nature of Applicant's Invention; Examiner's Cited References: 

As described in Applicant's originally filed application, in conventional hot deuterium 
fusion, the means for overcoming the Coulombic repulsion of two fusing nuclei is to provide one 
nucleus with a sufficient kinetic velocity with respect to the other. When the two nuclei are a 
critical distance apart, such distance being a function of their atomic number and atomic mass, 
the strong nuclear force dominates over the Coulombic repulsion and the nuclei fuse. 
Applicant's invention presents a method whereby a large surface charge on a conducting lattice 
effectively shields the positive charge of the one nucleus against the positive charge of the other 
so that a lower kinetic velocity is required to achieve the critical distance. Lower kinetic velocity 
translates to lower mean temperature. 

In cold fusion cells, the basic premise is that the unaltered intrinsic electron density of the 
metallic lattice provides this electronic shielding. With this limitation, cold fusion does not 
work, as is noted by Examiner, and as is referenced in the majority of papers provided by 
Examiner. However, one of the keys to Applicant's invention is a technique for significantly 
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enhancing the electronic shielding - a technology and methodology absent from all of the patents 
and papers referenced by Examiner. As an analogy, the late- 1800 car manufacturers did have an 
internal combustion engine. However, it was not until the weight-to-horsepower ratio of the 
engine was significantly decreased that the Wright brothers were able to use it to power their first 
airplane. That is, many fundamental inventions rely on known mechanical models that are 
optimized and better than those available in the past, crossing an operational threshold that 
makes the invention possible. Such is the case of Applicant's invention. 

Examiner references 14 patents and 107 non-patent documents. In general, the non- 
patent documents concern cold fusion, the basic premise of which Applicant already rejects . A 
year ago one may have made the argument that thermonuclear fusion "on a desktop" was a 
fantasy, but verified results from UCLA have refuted this argument (see Exhibit A). Clearly, 
fusion is slowly being removed from the fission thermonuclear bottle. 

As noted above, the quality of the shielding obtainable makes Applicant's invention 
possible. Such electron densities required can only be obtained through a charge concentrator 
with a precisely defined tip in the nanoscale. Just as the existence of macroscopic projections 
and protuberances on the electrodes are insufficient, so too will the Wright brothers' airplane not 
fly with an overweight and underpowered 1880 IC engine. 

Accordingly, recognizing the need for nanofabrication, one can dismiss all previous 
patents before about 1996, as such techniques were not available. The only patent referenced by 
Examiner after 1996 is U.S. Patent No. 6,489,704, which discloses a thermionic energy converter 
with fractional surface contacts, and whose technology and claims bear no relevance to the 
present invention. 
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Although Applicant has made his election of the claims, Applicant respectfully notes that 
Examiner has not recognized the common features in all of the groups of claims enumerated; that 
is, of a specialty surface generating sufficiently high charge density on the surface protuberances 
to allow significantly lower kinetic temperatures for hydrogen and its isotopes to nuclear fuse. 
This feature is described in the example given in Applicant's specification on page 23. The 
active portion or tips of the protuberances of Applicant's reacting surface are within tens of 
nanometers and, thus, are atomically sharp. 

Indeed, and contrary to Examiner's assertions, the manufacture of a nano-protuberance 
surface may be done, and has been done, by many known methods. In particular, scientists have 
constructed single protuberance structures to the single atom level, and have used them 
extensively in the tips of atomic scanning and probe microscopes. Applicant encloses herewith, 
for Examiner's reference, a report that Applicant has written describing such devices {see Exhibit 
B). Applicant further notes that there are companies commercially producing the atomically flat 
surface needed for fabrication - Applicant submits herewith literature from one such company 
(see Exhibit C). In further support of such nano-fabrication, Applicant submits an article 
showing a method for nanofabrication of a two dimensional array {see Exhibit D). Still further, 
there are numerous other methods in use for array fabrication {see Exhibit E, Fig. 2), including 
the techniques of nanolithography, and strain-induced pattern formation. Applicant encloses 
another article illustrating the growth of whiskers on a nanoscale {see Exhibit F). Applicant 
notes that the materials science literature in the last 7 years is replete with examples of such 
techniques of nanofabrication; these techniques are widely published and indeed known to the 
scientific community. In further support of such known technology, Applicant submits herewith 
a detailed engineering analysis completed to calculate the electric field and charge density on an 
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optimized cone separated from an electrically neutral flat plate, illustrating a practical 
construction rather than the optimized construction (zero separation distance) noted in the patent 
specification (see Exhibit G). 

On page 8 of Examiner's Office Action, Examiner refutes Applicant's total invention and 
takes no note of the theory presented in Applicant's specification. Examiner extensively quotes 
"accepted theory" without specifying how the snippet of theory presented varies from it. 
Applicant submits that the application is not intended to present a learned thesis on quantum 
mechanics. However, in the interest of explaining Applicant's Equations 1-3, Applicant encloses 
a small section on basic physics of the Bohr radius, indicating how the muon binding energy is 
calculated (see Exhibit H). From this expression for Bohr radius, clearly forcing the equality of a 
Bohr radius of increasing mass "m" and charge "Z" yields: 

0 0 0 0 0 0 

A z = h /(kZ e me) = A mu on = h /(ke m muon ) so that Z m e = m mu0 n, as noted. 

Applicant's equation considers an equivalent Z if all electrons were placed in a single 
energy shell. Although this violates the Pauli exclusion principle if the electrons are bound to 
the ion, it is quite different if the electrons are attached to the atomic metallic lattice on the 
surface of the cone, and the deuterium ion comes arbitrarily near this electron concentration. In 
this case, the solution of the Schroedinger equation for bound electrons on the deuterium ion 
nucleus does not apply. 

Examiner further objects to the form of Applicant's Equation 3 on multiple grounds, not 
the least of which is that the undetermined constants are not specified a priori. In response to 
Examiner's objection, Applicant submits herewith an interesting paper from The Astrophysics 

25 



Journal that looks at Thermonuclear Fusion and screening, and notes the fusion-enhancement 
factor due to screening in the solar plasma {see Exhibit I). 

Applicant notes that Equation 3 does not have a rigorous derivation, but stems from an 
observation that absent muon data, there is a reasonable straight-line correspondence between the 
logarithm of fusion temperature and the number of nucleon interactions for single proton fusion. 
The mathematical model is extended with the addition of the equivalent shielding created by the 
muon, and its corresponding fusion temperature. The coefficients are fit to the data by least 
squares, a rigorous mathematical method, which will yield only one solution (contrary to 
Examiner's assertions). All mathematical models of physical events must make assumptions to 
remain finite. The objective is to make the minimum number of assumptions so that the model 
may reasonably predict future behavior. 

Examiner further quotes Lindley who describes the electron cloud model in a metal with 
respect to elevation of mean electron mass. Such derivations of elevated electron masses more 
properly belong in BCS superconductor theory than in normal metals. No claim is made by 
Applicant that the electrons on the cone tip detach to the deuterium ions. Indeed, the entire 
concept of electron tunneling, used in numerous other inventions, relies upon the quantum 
mechanical uncertainty of the electron's exact position. Thus, many electrons may momentarily 
be arbitrarily close to the deuterium nuclei. Fukai is correct; deuterium nuclei in an unaltered 
palladium lattice cannot come close enough for fusion at low temperatures. However, this 
invention describes a highly enhanced surface charge density than that encountered in the 
unaltered metallic lattice alone. Like the Wright brother's IC engine, the enhanced charge 
density allows the fusion process. 
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Examiner further indicates that Applicant has provided inadequate disclosure of how 
Applicant's cones or wedges of the reacting surface can be manufactured - despite the fact that 
Applicant has already stated in the application that such manufacture may be had via 
nanotechnology. As discussed hereinabove, and as evidenced in the attached literature of the 
Exhibits, many known nanotechnology methods may be utilized to manufacture Applicant's 
invention. Applicant notes the recent textbooks of such methods, as referenced in Applicant's 
Exhibit B. Regarding electroplating, one of the primary mechanisms of failure of lead acid 
batteries is the growth of dendrites in the nanoscale between anode and cathode in the 
electrolyte. This growth has been studied for many years and its science is well developed. 
Further, there are multiple techniques in materials nanoscience including Atomic Layer 
Deposition (ALS) that can be utilized to provide precise and repeatable placement of atomic 
layers on surfaces. 

Applicant's Claim Amendments: 

To better clarify Applicant's invention, Applicant has now amended certain Independent 
Claims to reflect that each protuberance/cone/wedge of Applicant's reacting surface comprises 
an "atomically sharp apex." Applicant respectfully submits that the cited references and prior art 
do not teach or fairly suggest Applicant's invention as now claimed. In view of the foregoing 
arguments and evidence, reconsideration is respectfully requested. 
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CONCLUSION 



The above-made amendments are to form only and thus, , no new matter was added. 
Applicant respectfully believes the above-made amendments now place the Claims and 
application in condition for allowance. Should there be any questions or concerns, the Examiner 
is invited to telephone Applicant's undersigned attorney. 



Myers & Kaplan 
Intellectual Property Law, L.L.C. 
1899 Powers Ferry Road 
Suite 310 

Atlanta, GA 30339 
Phone: 770-541-7444 
Fax: 770-541-7448 



Respectfully submitted, 



Dated: September 





Ashish D. Patel 
Attorney for Applicant 
Reg. No. 50,177 
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Abstract 

Unique visualization of the true atomic surface structure is made possible by the 
geometric scanning and resolution of the scanning tunneling microscope (STM). In this device, an 
atomically sharp tip is maneuvered over the object of interest, with a potential being applied between it 
and the object. Instantaneous measurement of current in this circuit can be computer resolved into a 
simulated surface image. In practical design, the quantum electronic states of these tips are critical to their 
operation. This paper will review the literature on the theory of such tips in STM, and discuss 
considerations in their construction. A possible solution utilizing nanotechnological construction 
technique for problems encountered in their use will be presented. 

I. Introduction 

Unique visualization of the true atomic surface structure is made possible by geometric scanning and 
computer image reconstruction of the scanning tunneling microscope. The theory of scanning tunneling 
microscopy (STM) is reviewed quite well in Chapter 6 of the Woodruff and Delchar Text, Modern 
Techniques of Surface Science - Second Edition 1 . In brief, an atomically sharp tip is constructed and is 
controllably maneuvered over the object of interest, a potential being applied between the probe tip and 
the tip and the surface of the object. As a result of the electrical potential, a current is generated, the 
magnitude of which is a function of the separation distance from the tip to the surface and other factors. 
The exact position of the tip is controlled by the application of electrical potentials to mutually orthogonal 



t 

piezoelectric elements attached to the tip. The mechanical configuration is shown in Figure 1 , indicating 
the tip and piezoelements. 




Figure 1 . Basic Configuration of the Scanning Tunneling Microscope 
(taken from Figure 1 of Reference 2, Frantsesson and Zuev) 

In operation, the tip potential is kept low to prevent ion generation as is used in filed emission 
microscopy. The tip characteristics are critical, as the interaction of atomic orbital wave functions of the 
tip and the object surface determine the exact conductance. Thus, the exact material of the tip as well as 
its physical geometry with microscope operation must be considered. 



II. Quantum Theory Considerations 

The tunneling current, I, is proportional to the exponent 8 raised to the power -2KD, as in Equation II. 1 
where D is the effective distance between the tip and the object surface and K is the characteristic 
exponential inverse decay length. 



Eq.H.1 7 = C, -e~ 2KD where K = ^- 



where m is the electron mass, h is Planck's constant and <)> is the effective local work function. 1 

As there is some uncertainty in D and <)>, scientists have contemplated numerous models to remove the 
specific properties of the tip from the problem. In early work, Tersoff and Hamann considered the "s- 
wave model" in which the conductance is proportional to the local density of states p(r t ,E F ), as given in 
Equation H.2. 



Eq. H.2 p(r„E F ) = Y\<p v (r)\8(E v - E) r = r„E = E F 



where q> is the electronic wavefunction at the tip, E is the energy, and r is the position coordinate, taken at 
r t , the tip position. 3 In this model, the tip current maps out a contour of the constant values of the local 
density of states with the Energy E taken at the Fermi Energy level, E F . Unfortunately, this model yields 
resolutions no better than 6-9 Angstroms. 

A better model considers the P and D electron orbital states of the metal tips. This theory, advanced 
by Chen in 1990, yields substantially better resolution at a larger tip to sample distance. 4 However, to gain 
this resolution requires specific tip materials such as the D-band transition metals such as tungsten, 
platinum and iridium or certain semiconductor like silicon which form Pz-like metallic dangling bonds. 



III. Tip Fabrication 

By the Chen model, there are limitations as to the specific materials that can be used for the tips. 
Further, one must consider the method of manufacture of such tips. In the design of field ionization 
microscope tips one desires a nearly hemispherical and atomically smooth tip surface. One method of 
such construction is the application of a sufficiently high electric field to the rough tip, which blasts any 
protrusions from the tip as positive emitted ions. Solution of Maxwell's equations in the static regime in 
the region around the sharply pointed tip indicates a major charge concentration effect. As nature tends 
always to smooth out irregularities, then a smooth hemispherical surface will be created. However, this 
result is not appropriate for the resolution desired in the STM. 

Aeshimann et al describes a process of creation of silicon oxide tips, as shown in Figure 2. These tips, 
utilized for scanning near-field optical microscopy utilize 12-micron silicon oxide tips integrated over the 
end of a silicon cantilever. 5 Essential to this design is repetitive micro machining to yield tip geometry. 




Figure 2. Silicon-dioxide tips (taken from Figure 1 of Reference 5, Aeshimann et al) 



As noted by Liu and Gamble, desirable properties for the scanning probe microscope tips include 1) 
good tip sharpness with small tip radius 2) small force constant and 3) high resonant frequency. 6 In their 
process, the sharp tips are formed by oxidation sharpening rather than direct anisotropic etching (AE). In 
AE, a photolithographically patterned region of silicon dioxide forms on an epitaxial silicon layer that is 
then immersed in ethylene diamine pyrocathechol (EDP) producing pyramid shaped tips. The authors note 
that while this process, AE, may produce atomically sharp tips, it is difficult to control due to variations in 
process parameters. The authors improve the general process by serial growth of thermal oxide on the tip 
unit and subsequent removal, in a non-uniform way, of the oxide using dilute hydrofluoric acid. By this 
subtractive process, the apex becomes sharper after each sharpening cycle. 

A favored technique for tip fabrication involves plasma etching and oxidation sharpening, as reviewed 
by Tabib-Azar and Wang, and is shown in Figure 3. In their process, the starting substrate is a double- 
side-polished silicon on insulator (SOI) wafer with a 1 5-uM device layer, a 1-uM buried oxide layer and a 
400-uM handle layer. After cleaning the wafer has a thermally oxidized 1-uM silicon dioxide layer grown 
which is photolithographically fashioned into defined circular patterns. The exposed thermal oxide is 
etched by buffered hydrofluoric acid and then silicon fluoride plasma etchings produces conical silicon 
tips, as shown in Figure 3 A. Conventional thermal oxidation sharpening produces a peak as shown in 
Figure 3B, with a 50 Angstrom radius. Finally, in their process, a coaxial shielding layer is built up 
around the tip, as shown in Figure 3C, for better STM electrical performance. 7 
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Figure 3. Tip SEMS (taken from Figure 8 of Reference 7, Tabib-Azar and Wang) 



In a more advanced fabrication methodology, Cruchon-Dupeyrat et al utilizes CAD design techniques 
and automated vector scanning to do the nanofabrication of the tip 8 . In their process, a multi-mode 
scanning probe microscope (SPM) is linked with a computer control for precise automated control of the 
SPM tip. This arrangement allows almost molecule-by-molecule fabrication of the new tip or indeed any 
desired nanostructures. While offering unprecedented control, it suffers from the constant defect of all 
devices of this class, that is, the excessive time required for complete fabrication before contamination 
sets in, even under ultra high vacuum conditions. It would be preferable to use nanotechnology self- 
assembly methods rather than atom be atom placement. 



IV. Problems in Current Production Methods 



At present, probe tips are constructed on a one by one basis. Built up over a conventional substrate, 
they experience constant contamination and, with loss of sharpness and concurrent degradation of STM 
performance, a need for replacement. Further, there is difficulty in cleaning the tips, a process to remove 
secondary adhered species without losing the atomic sharpness. Conventional fabrication is a subtractive 
process, with etching cycles. Simply applying a high potential to "blast-off 9 secondary adherents would 
work only if the exact tip geometry were known and uniform so that the applied potential could be kept 
low enough to prohibit loss of the metallic tip atoms themselves. Unfortunately, variation in the 
parameters of manufacture prevents this cleaning approach, with the result that frequent complete tip 
replacement is required for proper STM operation, leading to high operating costs. Further, each tip 
replacement may require instrument recalibration due to non-homogeneity of the tips. 

V. Proposed Nanotechnological Engineering of Tips 

Fundamental to solution of uniformity in tip construction is uniformity of construction of the tip 
substrate. Known nanotechnological materials processing techniques are used to create an atomically flat 
substrate, and the tip can be created in an additive fashion rather than a subtractive fashion, allowing easy 
removal of adhered impurities and renewal to atomic sharpness. 

An atomically flat substrate can be created by essentially "freezing" a material in a non-equilibrium 
state. If a molten metal is cooled at a rate exceeding 1 million degrees Centigrade per second, then normal 
formation of material dislocations and voids are severely inhibited, with the result that an atomically flat 
state results. Currently in commercial use, the chill-block-melt-spinning apparatus is used to create an 
atomically flat metal ribbon. Powdered metal is injected under pressure into a cylinder and the metal is 
then melted. A nozzle from the cylinder is directed upon a rapidly spinning cooled wheel and an atomized 
metal is then sprayed on this moving surface with the result that a thin ribbon of metal is produced 
substantially free of surface roughness and defects. 



From the flat substrate one may then build up a conical microstructure. One method is use to create 
multiple peaks, whose number may be regulated by surface masking is by laser interference metallic 
deposition. As explained by G. Cau in Nanostructures and Nanomaterials a nd shown in Figure 4, two 
laser beams are placed parallel to an atomicallyjlat surface and perpendicular to each other, creating an 
electric field interference pattern. 9 Placement of this field enhanced surface in a vapor of metal ions will 
cause such ions to deposit and align according to the interference pattern leading to a highly uniform and 
regular multiple tip unit. With uniform geometry and materials, the variations in manufacture are 
minimized, and with multiple close tips instead of one creating the STM signal, D and P orbital patterns 
uniquenesses of single peaks are minimized. 



oven 




optical potential 




h-X/2 -< 



distribution of 
atoms 
ontbe substrate 

(simulation) 




Schematic illustrating the basic principles of neutral atom lithography with light 
forces. (B. Brezger, Th. Schulze, U Drodofsky, J. Stuhler, S. Nowak, T. Pfau, and 
IMiynek, 1 Vac, Sci Technoi B15, 2905 (1997).] 



Figure 4. Neutral Atom Lithography (taken from Fig.7.7 of Reference 9, Cao) 



Other methods of construction of micro cones on atomically flat surfaces can be used. For example, one 
could laser drill an atomically flat mask, place such a mask over an atomically flat substrate and then 
create regularly placed growth nidus by vapor deposition of metal atoms in a clean atmosphere on the 



substrate through the holes. Placement of this modified planar surface in an electrolytic bath, with a 
second non-modified atomically flat surface facing it, and applying an electric current with the correct 
intensity and polarity will cause electrolytic deposition of atoms of the second surface preferentially upon 
the growth niduses, the locations of least electrical resistance. In this way a set of atomically sharp conal 
whiskers are created, in an array or arbitrary size, with uniform properties for the composite STM tip. 
Further, in the event that these multiple cones become contaminated during STM operation, one could use 
high voltage blow-off of the contaminant, as the underlying tip structure is well characterized. Even if 
there were some loss of sharpness in the cleaning process, one would reinsert the blunted composite tip 
into an electrostatic bath as noted to regrow the atomically sharp conal whiskers. Use of a 2x2 or 3x3 
multiple tip matrix would also increase STM sensitivity by providing a scanned signal as the object 
surface protrusions are brought under sections of the tip peaks. Further, as opposed to using piezoelectric 
manipulation to control probe tip movement in the orthogonal x,y,z planes only, one could place such 
piezocrystals at oblique angles to the surface allowing angular control of the composite tip, and tracing 
over the descending slope of the object surface hillocks. This control may be though of as similar to the 
head motion of a 6-axis CNC drilling machine as opposed to the existent x,y,z 3-axis movement. 

VI. Conclusion 

Current methods of tip manufacture, primarily subtractive, while creating atomic sharpness necessary 
for STM operation, suffer from non-uniformity and lack of ability to be effectively cleaned. Constant 
replacement is costly, both for the actual replacement tip, and for the need for STM recalibration with 
each new tip. Further, single sharp tips for STM require some empirical instrument calibration as the 
exact geometry of the tip is not known a priori. The proposed new method of manufacture is an additive 
process, allowing highly uniform tip fabrication, manufacture of multiple uniform tips on a single 
substrate, allowing multi-tip scanning to minimize the geometry effects of a single tip, and an additive 
process for tip reconstitution after elimination of impurity adherents. 
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dense parts in a single process in any form factor required. Sheet-stock of our toughest 
materials are now available. 

Note: All prices listed in US Dollars 




Other Materials: 

* Nano Coatings 

♦ Nano Powders 
GBE® Materials 



Manage Order 



Nickel (99.9%) 

Attributes: -20nm average grain 
corrosion resistance. 
Applications: Springs, MEMS, high 

ttem# 



NS-F-1A 
Ni-F-1 B 
Ni-F-2A 
Ni-F-2B 
Ni-F-3A 
Ni-F-3B 
Ni-F-4A 
Ni-F-4B 
Ni-F-5A 



Sample Size 
in (mm) 

0.002 x 1 x 1 
(0.050 x 25.4 x 25.4) 

0.002 x 3 x 3 
(0.050x76.2x76.2) 

0.004x1 x1 
(0.10x25.4x25.4) 

0.004x3x3 
(0.10x76.2x76.2) 

0.006x1 x1 
(0.15 x 25.4 x 25.4) 

0.006 x 3 x 3 
(0.15x76.2x76.2) 

0.008x1 x1 
(0.20x25.4x25.4) 

0.008x3x3 
(0.20 x 76.2 x 76.2) 

0.02 x 1 x 1 
(0.50x25.4x25.4) 



- high strength, high resilience, excellent 

yield sputtering targets 

Unit Price Quantity 
Required 



Orde r 



$115.00 
$750.00 
$120.00 
$800.00 
$130.00 
$850.00 
$135.00 
$900.00 
$300.00 
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1 
1 
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Nickel- 20%Iron 

Attributes: -20 nm average grain size - high strength, high resilience, excellent 



magnetic properties. 

Applications: Springs, electrical contacts, s< 


rft magnets 








ltem# 


Sample Size Unit Price 
in (mm) 


Quantity 
Required 


Manage 
Border 


Ni20Fe-F-1A 


0.002 x 1 x 1 
(0.050x25.4x25.4) 


$130.00 


1 




Add 












Ni20Fe-F-1B 


0.002 x 3 x3 
(0.050 x 76.2 x 76.2) 


$850.00 


1 




Add 


Ni20Fe-F-2A 


0.004x1 x1 
(0.10x25.4x25.4) 


$135.00 


1 




Add 


Ni20Fe-F-2B 


0.004x3x3 
(0.10 x 76.2 x 762) 


$900.00 


1 




| Add 


Ni20Fe-F-3A 


0.006x1 x1 


$150.00 
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Decreasing grain size, increasing performance; Integrals 
Bottom-Up Nanomaterial Process delivers unparalleled material 
performance in a cost effective one-step process. 

Via a patent-protected electroeynthesis technology, materials having average crystal sizes 
1000-Bmes smaller than those found in conventional materials can be inexpensively 
produced. 

Decreasing the Grain Size of Materials... 

Tvoical metals and ceramics have crystalline structures with average grain sizes of several 
microns. Integrals patented bottom-up process, creates materials with 1000 times smaller 

NNHIMIBMl PROCESS 1 ^^^.^^^^.^ 

Conventional Material Nanocrystalline Material TEM Magnified view 

p^i^iQtirAfioilcattons 

Order Manomaterials Now 
• Powders 
- Coatings 
Polls and Plates 
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The nanocrystalOne material has an average grain size of 20 nm. This is about 100Qx 
smaller than the conventional material, at about 20 microns 

...Drastically Improves Material Performance 

Properties of metals are governed by the Hal l-Petchi ^°^'^J^^£^^ 
strength increases, hlanocrystaltine materials are characterized with significant increases 

in: 

■ Yield Strength 

■ Ultimate Tensile Strength 

■ Hardness , 
Without Sacrificing Ductility t , ... 

: A value of 5% elongation; often used by structural designers, is used to separate ductile 
from brittle behavior. To date. (Integrals) nanomaterial may be the ^jy exampte of 
nanocrystalline material produced by a one-step process exhibiting ductility behaviour 

C.C. Koch in "Ductility in Nanostructured & Uttrafine Grained Materials: Recent Evidence for Optimism" 
J. Metastable & Nanocrystalline Materials. Vorf. 18 (200$pp.9-20. 

Other characteristics of nanocrystaffine materials include: 

■ Wear Resistance & Resilience 

■ Thermal Shock and Fatigue Resistance 

■ High Magnetic Saturation and Coercivity (@60Hz) 

■ Enhanced Anti-Microbial Activity 



http://wwwintegranxom/tech/process.hto 



2/11/2005 



PATENT 
10/735,406 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 



In re: Robert Indech 

Serial No.: 10/735,406 

Filed: December 12, 2003 



Docket No: 23050-RA 
Examiner: Behrend, Harvey E. 
Group Art: 3641 



For: APPARATUS AND METHOD FOR FACILITATING NUCLEAR FUSION 

RESPONSE TO OFFICE ACTION OF JULY 19. 2005 

EXHIBIT D 



Nanofabrication of a two-dimensional array using laser-focused atomic 
deposition 

R. Gupta, J. J. McClelland, 8 * Z. J. Jabbour, and R. J. Celotta 

Electron Physics Group, National Institute of Standards and Technology, Gaithersburg, Maryland 20899 
(Received 16 May 1994; accepted for publication 27 June 1995) 

Fabrication of a two-dimensional array of nanometer-scale chromium features on a silicon substrate 
by laser-focused atomic deposition is described. Features 13± 1 nm high and having a full-width at 
half maximum of 80± 10 nm are fabricated in a square array with lattice constant 212.78 nm, 
determined by the laser wavelength. The array covers an area of approximately 
100 /xmX200 /xra. Issues associated with laser-focusing of atoms in a two-dimensional standing 
wave are discussed, and potential applications and improvements of the process are mentioned. 



New methods for fabrication of nanometer-scale struc- 
tures have been under intensive investigation recently be- 
cause of the perceived benefits which might arise, such as 
smaller electronic devices, higher-density information stor- 
age, and novel materials. Within the past three years, a new 
technique for nanostructure fabrication involving laser- 
focused atomic deposition has been demonstrated. 1 " 3 In this 
paper we present a significant enhancement of the original 
technique, which was used to fabricate lines on a substrate, 
to demonstrate the first two-dimensional fabrication. With 
these new results, we also discuss some of the considerations 
associated with generalizing the dimensionality of the pro- 
cess. 

In laser-focused atomic deposition, a laser light-field is 
used to control the motion of atoms as they deposit onto a 
surface. This approach has a number of potential advantages 
in comparison with conventional fabrication techniques, such 
as optical or electron-beam lithography. The advantage over 
optical lithography lies in the potential for much higher reso- 
lution. Optical techniques are fundamentally limited by dif- 
fraction to a minimum feature size of about half the wave- 
length of the light used. For ultraviolet light, this corresponds 
to about 100 nm. Free-flying thermal atoms, on the other 
hand, have De Broglie wavelengths of order 10 pm, so dif- 
fraction effects can in principle be reduced to a negligible 
level. The actual resolution limits of laser-focused atomic 
deposition are still a subject of research, though structures 
have already been demonstrated at the 65 nm level and theo- 
retical predictions suggest that 5-10 nm features may be 
possible. 2 

This range of feature size is already just attainable with 
electron beam lithography; 4 however this process is inher- 
ently serial, that is, a complex pattern must be fabricated by 
scanning the electron beam across the surface. For large, 
complex patterns, the fabrication time and associated drift 
problems make electron beam lithography less desirable. 
Laser-focused atomic deposition, on the other hand, does not 
suffer from these limitations because it can be implemented 
in a parallel fashion. 

A one-dimensional schematic of the two-dimensional 
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laser-focused atomic deposition process used in the present 
work is shown in Fig. 1. A laser standing wave is generated 
across the surface of the substrate, and chromium atoms, 
collimated to 0.25 mrad in each of two dimensions by laser 
cooling, 5 are directed at the surface, traversing the laser field 
on their way to deposition. The laser field, produced by a 
single-frequency CW dye laser, is tuned 500 MHz (100 natu- 
ral line widths) above the atomic resonance line in chromium 
at \ = 425.55 nm (vacuum wavelength). With this tuning, a 
dipole force 6 is exerted on the atoms toward the low intensity 
regions of the light field. The result is a concentration of 
atoms at the nodes of the standing wave, which occur at 
intervals of X/2= 212.78 nm. 

In one dimension, the process as depicted in Fig. 1 is 
relatively straightforward. One aspect of the geometry illus- 
trated in Fig. 1 that simplifies the implementation of this 
process is the fact that the positions of the deposited lines on 
the substrate depend, only on the standing wave node posi- 
tions, which in turn depend primarily on the absolute dis- 
tance along the substrate from the mirror generating the 
standing wave. To first order, variations in the position or 
direction of the laser beam have no effect on the node posi- 
tions, and as long as good stability is maintained between 
mirror and substrate, the periodicity of the pattern will be 
stable. 



Chromium atoms 

wuwwwwwwwhuhw 




213 nm 



FIG. 1. One-dimensional schematic of laser-focused atomic deposition pro- 
cess, showing chromium atoms being focused by a laser standing wave into 
its nodes. The trajectories and the deposited peaks represent the results of 
actual calculations of the focusing process, though the relative vertical 
scales are highly distorted for clarity. 
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For a two-dimensional deposition, the standing wave is 
formed by crossing two one-dimensional standing waves at 
90° across the substrate. With this higher dimensionality, the 
nodes of the two-dimensional standing wave become more 
complex, depending on the polarizations of the constituent 
standing waves and also in some cases on their relative tem- 
poral phase. There are still zero points in the light intensity at 
locations on the surface that are an integral number of half- 
wavelengths from each standing wave mirror, but examina- 
tion of an expression for the net electric field amplitude 7 
shows that additional nodal patterns can exist, and these can 
shift around as the temporal phase difference between the 
standing waves varies. 

One approach to eliminating unwanted variations in the 
nodal pattern due to temporal phase variations is to generate 
the two orthogonal standing waves in an actively-stabilized 
optical cavity. 8 Another approach, 9 which is useful when sta- 
bility of the nodal topography is desired but absolute position 
is unimportant, makes use of three laser beams intersecting 
at 120°. For the present work we have made use of the fact 
that one can choose a polarization scheme in which the tem- 
poral phase does not affect the nodal pattern or position. If 
one chooses orthogonal linearly-polarized standing waves, 
with one polarized perpendicular to, and the other parallel to, 
the plane of the substrate, the resulting intensity distribution 
is given by 

I(x,y) = 4/ 0 ( sin 2 kx + sin 2 ky) ( 1 ) 

where I Q is the intensity of a single incident wave, 
k=2ir/\ is the wave vector of the light, and the x- and 
y-axes are defined to lie in the plane of the substrate. The 
intensity given by Eq. (1) forms a pattern of nodes and peaks 
on a X/2XX/2 square lattice with the peaks separated by 
saddle regions along i and y at half the intensity. 

While the intensity of Eq. (1) does not suffer from tem- 
poral phase disturbances, it has potential drawbacks that 
could in principle cause problems with a deposition experi- 
ment, though these appear to be less severe in practice. First, 
it might seem that since the intensity does not have cylindri- 
cal symmetry around the potential minima, the resulting de- 
posited spots would not be round. However, the intensity is, 
in fact, surprisingly symmetric in the regions near the 
minima. This can be seen mathematically by converting x 
and y to the polar coordinates (r,0) and noting that 
I(r,6)**4I 0 k 2 r 2 (i.e., the ^-dependence drops out) for 
kr<\. Exactly how much effect the non-symmetric regions 
of the potential (where kr is not much less than 1) have on 
the deposited pattern relative to the symmetric regions is 
difficult to predict without a ray-tracing calculation. 10 Nev- 
ertheless, the experimental results indicate that there is little 
effect. 

Second, it must be noted that while the intensity in this 
configuration does not depend on the relative temporal 
phase, the polarization of the electromagnetic field is com- 
plicated. In fact, the local polarization varies dramatically as 
a function of x and y over the scale of a wavelength, and 
exactly what form this variation takes depends on the relative 
temporal phase. The local polarization of the field can be 
very important because it can determine the strength of the 
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FIG. 2. Atomic force microscope (AFM) image of chromium features 
formed by laser-focused atomic deposition in a two-dimensional standing 
wave. The image covers a 4/imX4/im region of the sample. The features 
are on a square lattice with spacing 212.78 nm, which is determined by the 
laser wavelength. The standing wave is formed by superimposing two one- 
dimensional standing waves oriented at 45° and 135° to the figure. Also 
shown are two line scans, labeled A-A' and B-B', whose locations are 
indicated on the AFM image. The vertical scale for the line scans was 
determined by the AFM calibration, and the offset was estimated by requir- 
ing that the integral under the surface equal the total amount of material 
deposited, obtained from flux measurements. 

laser-atom interaction. Chromium atoms entering the 
standing-wave field are in their ground state, but are distrib- 
uted among a number of degenerate magnetic sublevels. The 
strength of the laser-atom interaction for each magnetic sub- 
level varies for different polarizations of the laser, due to 
differences in the Clebsch-Gordan coefficients for the vari- 
ous transitions between the different magnetic sublevels. 11 
The result is a potentially wide variation in the force on the 
atom as a function of space, and also time if the temporal 
phase is not stabilized. 

Despite these potential problems, it appears that a suffi- 
ciently symmetric and stable potential exists for a well- 
defined pattern of features to be created. Figure 2 shows an 
atomic force microscope (AFM) image of a 4 /imX4 /xm 
region of the two-dimensional chromium pattern deposited 
on a silicon substrate at room temperature. For this deposi- 
tion, the laser beams for the two dimensions had Me 1 diam- 
eters of 0.13 ±0.02 mm and each contained a single-beam 
traveling-wave power of 12± 1 mW. 12 The total deposition 
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time was 20 min. The full pattern covers an area of about 
100 /imX200 fim, and is quite uniform across the entire 
area. The lattice constant of the pattern, \/2 or 212.78 nm, 
was considered to be known more accurately than the hori- 
zontal calibration of the AFM, so this value was used to put 
the horizontal dimensions of the image on an absolute scale. 

Also shown in Fig. 2 are two line scans, showing the 
shape of the features in two directions. These line scans have 
been put on a true vertical scale, including the background 
level, by normalizing the surface topography to the average 
deposition thickness of 20 nm, estimated (with an accuracy 
of about ±5 nm) from previous characterizations of the 
deposition rate. Line scans such as the ones depicted in Fig. 
2 and others taken at a variety of locations on the sample 
indicate that (without correction for AFM tip shape) the fea- 
tures are 13± 1 nm high, and have a full-width at half maxi- 
mum of 80± 10 nm. The cause of the slight asymmetry in 
line scan A-A' is unknown, though we believe it to be an 
AFM artifact. 

As the line scans indicate, the regions between the fea- 
tures are also covered with chromium. This background is a 
result of other isotopes in the atomic beam which are not 
affected by the laser light (16% of the atoms), atoms which 
are transferred to the metastable D-state during the optical 
collimation step (an estimated 7% of the atoms), and atoms 
that are in the high velocity tail of the Maxwell-Boltzmann 
velocity distribution emerging from the chromium evapora- 
tion source. Scan B-B' appears to have a slightly higher 
background level, arising we believe because it samples the 
saddle regions of the standing wave, while A-A' samples the 
true nodes. 

Various mechanisms could be implemented to reduce the 
background seen in Fig. 2, if desired. For example, the un- 
desirable isotopes and D-level atoms, as well as some of the 
high velocity atoms, could be removed from the beam by a 
laser deflection process. Alternatively, it is possible that a 
uniform etch of the chromium surface as deposited could 
remove the background before eliminating the features com- 
pletely. 

This work demonstrates that laser-focused atomic depo- 
sition can be successfully used to create a uniform, two- 
dimensional nanometer-scale pattern on a substrate. In its 
current form, the pattern shown in Fig. 2 could prove ex- 
tremely useful as a calibration standard on the nanometer 
scale. The features represent essentially a "contact print" of 
a light wave, the wavelength of which is tuned with ex- 
tremely high precision to an atomic resonance whose fre- 
quency is know with very high accuracy (about 1 ppm). 
Thus, within the limits of some geometrical corrections that 
can be kept very small (of order 10 ppm or less), the lattice 
spacing accuracy is extremely high. 



With additional enhancements of this technique, further 
applications can be envisioned. Elimination or removal of 
the background will result in the creation of an array of iso- 
lated nanoscale metal dots on a surface. These could be used 
to study, e.g., transport phenomena, or quantum dot effects if 
fabricated on a semiconductor. Further, they could be used as 
an etch mask 13 to transfer the pattern to a substrate material, 
allowing extension of the fabrication techniques to other ma- 
terials. Improvement of the resolution, which should be pos- 
sible down to the 10 nm level, could lead to the possibility of 
scanning the substrate during deposition, creating almost any 
desired pattern, replicated within each unit cell of the lattice 
across the substrate. Furthermore, the interference of many 
standing waves incident from a range of angles with con- 
trolled phase could be used to generate more complex pat- 
terns. 

This work is supported in part by the Technology Ad- 
ministration of the U. S. Department of Commerce, and by 
the National Science Foundation under Grant no. PHY- 
9312572. Z. J. J. acknowledges support of a National Re- 
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Field emission from silicon spikes with diamond coatings 
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Diamond-coated silicon field emitters were fabricated and investigated. Emission currents of few 
M per tip at voltages of several hundred volts were obtained from very blunt tips with curvature 
radii up to 3 /tm. The values of the effective work function calculated from Fowler-Nordheim plot 
were in the range from 0.3 to 1.2 eV. Two models for an explanation of the experimental data are 
proposed. © 1995 American Vacuum Society. 



I. INTRODUCTION 

Recent attempts in the application of diamond-coated tips 
as field emission sources demonstrated new emission phe- 
nomena that could be considered as very promising for 
vacuum microelectronic devices. There are several important 
advantages of diamond from the field emission point of view. 
(1) The negative electron affinity (NEA) of the (111) plane of 
diamond 1 has been suggested as the cause of observed low 
effective work functions obtained from a simple implemen- 
tation of the Fowler-Nordheim equation. (2) The thermal 
conductivity of diamond is 20 W/cm °C, that is, the highest 
of all known materials. 2 This property is favorable for high- 
current applications. (3) The excellent chemical stability of 
the surface of diamond could enable current stability and low 
vacuum operation. (4) High breakdown threshold (10 7 V/cm) 
and rather large high-field electron velocity (2.7 X10 7 cm/s) 2 
should enhance the maximum emission current, limited by 
current saturation. 

The first report about growth of polycrystalline diamond 
particles on the very end of silicon whiskers by a hot fila- 
ment chemical vapor deposition (CVD) method was pre- 
sented in Ref. 3. Bias enhanced microwave plasma CVD also 
was used for the coating of silicon tips. 4 In Refs. 4 and 5, 
preliminary data of field emission from "diamond" coated 
silicon whiskers were presented 4,5 shortly thereafter. Emis- 
sion phenomena such as high emission currents from very 
blunt tips, 5 large emission area, and high current stability 4 
were reported. 

Prospects are good for the application of diamond-coated 
tips for field emission devices. In this article new experimen- 
tal data and results of calculation of emission parameters 
using various approaches are presented. 

II. EXPERIMENT 

A. Preparation of the tips and diamond growth 

The initial tips were vapor-liquid-solid grown silicon 
whiskers. As-grown whiskers were of about 100 jjltu in 
height and several /urn in diameter. Each whisker had a 
"cap" of a Au-Si eutectic alloy, a few fim in size on the 
very end (Fig. 1). The tips can be sharpened by wet etching 
and oxidation techniques and the cap is removed by sharp- 
ening (Fig. 2). 

Diamond particles were grown on regular arrays of the 
tips by a hot filament CVD method. 3 Electron diffraction 
studies had shown that the particles consisted of polycrystal- 



line diamond with the average size of the crystallites about 
10 nm. 5 No investigations concerning band structure, de- 
fects, surface, and gap states were performed. Both as-grown 
samples with gold caps and sharpened ones were used in 
these experiments. The preferential growth of diamond at the 
very end of the tips, especially in the case of sharp whiskers, 
took place in these experiments. The possible explanation of 
this phenomenon was proposed in an earlier paper. 3 A regular 
array of silicon tips with diamond particles is shown in Fig. 
3. 

Depending on the preparation of the samples and the 
growth conditions, three kinds of tips were obtained: (1) 
blunt tips with a gold cap coated by a lot of diamond par- 
ticles (Fig. 4); (2) sharp tips (radii of curvature from 10 to 
100 nm) that had only the spherical diamond particle on the 
very end ("ball-on-point") (Fig. 5); and (3) sharp tips coated 
by a nearly continuous layer of coalescent particles of dia- 
mond ("coated point") (Fig. 6). 

The samples were prepared for field emission measure- 
ments by "shaving," i.e., by carefully removing (breaking) 
whiskers off the substrate, leaving the only tip on each 
sample in order to perform experiments under well-defined 
conditions. 

B. Field emission 

The geometry used for the field emission experiments is 
shown in the Fig. 7. The measurements were made in an ion 
pumped vacuum chamber at pressures of 10" 7 -10" 9 Torr. 
The anode was a stainless steel disk 50 mm in diameter. The 




Fig. 1. As-grown silicon whiskers with caps of Si-Au eutectic alloy. 
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Fig. 2. Sharpened silicon whiskers. 

spacing between the base of the sample and the anode of 
about double the height of the tip was maintained (see Table 
I). 

Three kinds of field emitters were investigated (see Sec. 
II A, Figs. 4, 5, and 6). Current- voltage (/- V) and Fowler- 
Nordheim (FN) plots for each of the three emitters are shown 
in Figs. 8(a) and 8(b). 

From these experiments high emission currents (a few fiA 
per tip) at rather low voltages (hundreds of V) were achieved 
from very blunt tips (up to several /-an). These results were 
obtained on a number of samples so they can be considered 
as consistent among similar tips. All emission data were 
summarized and published elsewhere. 5 

No pressure dependence of the emission current was ob- 
served in the range of 10" 7 -10~ 9 Torr. (A pressure of 10" 9 
Ton* was achieved first by baking in the chamber; leakage of 
air was used to change the pressure. No analysis of partial 
pressures in the chamber was performed.) 

III. RESULTS AND DISCUSSION 

The emission parameters were obtained from FN plots by 
the slope-and-intersect method. The conventional field emis- 
sion is described by the Fowler-Nordheim equation: 

J=(AE 2 /(<f> t 2 (y)))ex V [-B<f> 3/2 v(y)/El (1) 

where J is the current density, E the field strength at the tip, 
/the work function, ,4 = 1. 54 Xl(T 6 , £=6.87X10 7 , the func- 




Fig. 4. The very tip of a not sharpened Si whisker with a gold cap coated by 
diamond particles. 



tions t(y) and v(y) are usually approximated by t(y) = \.l 
and v(y)=0.95-y. By substituting I=aJ and E=/3V 
(where a is the emitting area and @ the field enhancement 
factor), Eq. (1) can be rewritten in the form obtained in the 
experiment: 



I=aV 2 exp(-b/V). 



(2) 



We used the point-to-plane geometry expression 6 for /? in 
our case (see Fig. 7): 



P = (h/r)(V/d). 



(3) 



The tip's height h and curvature radius r were measured 
by scanning electron microscope (SEM) with good accuracy 
thus abling us to calculate /3 directly. 

The values of the work function <f> and effective linear 
dimension of the emission area a, found from FN plots, are 
presented in the Table I. The.se data are just the result of 
formal application of the Fowler-Nordheim equation devel- 
oped for metals and for very approximate values of t{y) and 
v(y). No physical meaning could be extracted from this ap- 
proach. 

It should be noted that the FN constant A = 1 .54X 1 0~ 6 is 
derived for degenerated electron gas. In general, semicon- 
ductors may have a different value of this constant, depend- 
ing on carrier concentration. Since real concentration in elec- 
trons in our case is not known and the emission area, found 
from the FN plot depends on A, the formal calculation will 






Fig. 3. Regular array of silicon whiskers with diamond particles. Fig. 5. Sharpened Si tips with only the diamond particle on the very ends. 



JVST B - Microelectronics and Nanometer Structures 



420 Zhlrnov, Givargizov, and Piekhanov: Field emission 




3348 29KV X9, 600 I ^ ro WD 1 4 



Fig. 6. Sharpened Si tips coated with a nearly continuous layer of coalescent 
diamond particles. 



not give a real value of the emission area. However, there are 
several reasons why we used a "metallic" value of A . 

(1) Since the diamond coatings in our experiments be- 
haved like conductors and the saturation limit was not 
achieved, we assumed that the Fermi level position is near 
one of edges of the band gap. Further, since current satura- 
tion had not been reached in experiments, the "zero-current" 
approximation should be valid, i.e., the number of emitting 
electrons is considerably less than the number of electrons 
near the surface of the tip. A constant A is therefore assumed 
to be close its metallic value in this case. 

(2) The value of A = 1.54X10~ 6 is often used for w-type 
Si (see, for example, Ref. 7). A comparison of emission areas 
before and after diamond deposition can provide useful in- 
formation about the modification of emitter properties even 
if the change of "effective emission area" is due to a change 
of electron concentration. Thus, by the term "effective emis- 
sion area" we mean not the real emitting part of the tip but 
an additional parameter for the characterization of emission 
behavior. A comparison of noncoated and coated Si emitters 
for A = 1.54X 10" 6 has been done 4 previously and the authors 
reported a considerable increase in the effective emission 
area after diamond coating. 

In order to obtain more accurate data a "modified for a 
semiconductor case" FN equation was used. First, the ex- 




Fig. 7. Geometry of the field emission experiments (values d, h, and r are 
given in Table I). 
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Table I. Experimental emission parameters of diamond-coated field emit- 
ters. 



FN 



Emitter (fim) (jim) (jim) a b (cm) (eV) 

Blunt tip with 200 400 3 2.8X10" 9 8688 2.5XlO~ 8 0.37 
coating (Fig. 4) 

"Ball-on-point" 100 200 0.4 3.1XI0" 9 5372 7.9X10* 8 1.01 
(Fig. 5) 

"Coated point" 200 400 0.4 3.8X10' 9 4146 5.7X10 -8 0.86 
(Fig. 6) 



pression for Schottky lowering of the work function barrier 
y(f) contained a dielectric constant of material K in this 
case: 

y = 3J9X\0- 4 (E m /<f>)(K s _ x )/(K s + x ). (4) 

Second, the approximation of functions t(y) and v(y) in 
form of polynomials 7 was used. 

Two different hypotheses were under consideration: (1) 
real lowering of the effective work function; (2) observed 
emission behavior not due to a low work function of blunt 
tips but due to the rough morphology of the tips, i.e., emis- 
sion occurs from many of the very sharp microtips on the 
surface of the blunt tip. Then we consider both hypotheses in 
more detail. 

(1) Real lowering of the work function. Assum- 




Fig. 8. (a) l-V and (b) FN plot for three kinds of emitters: not sharpened 
whisker with gold cap coated with diamond particles (1), sharpened tip with 
only the diamond particle on the very tip (2), and sharpened tip coated with 
a layer of coalescent diamond particles (3). In (a) model curves, calculated 
using LWF and TSFE approaches are also given. 
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Table II. Emission parameters of diamond coated field emitters calculated 
using LWF and TSFE approaches. 





LWF approach 




TSFE approach 








h x r, 0 


*i 


r, h 2 r 2 






Emitter 


(/mi) (nm) (eV) 




(nm) (nm) (nm) 


N 


(eV) 


1 


200 3000 0.37 


200 


3000 12 0.2 


200 


5.4 


2 


100 400 1.03 


100 


400 8.7 0.7 


10 


5.4 


3 


200 400 0.88 


200 


400 9.5 0.6 


10 


5.4 



ing the real lowering of the effective work function (low 
work function (LWF) model), e.g., tunneling of the electrons 
from the conduction band or from the surface states near the 
conduction band, we calculated the I-V curves for A^=5.6 
which is known for diamond and for the geometry of the 
emitters used in our experiments (see Table II). 

(2) Emission from small protrusions on the surface. Cal- 
culations of the emission current for this case were made on 
the assumption that electrons that tunnel from the valence 
band mainly contribute to the emission current. The value of 
the work function, approximately equal to the band gap of 
diamond, 0=5.4 eV, was used in the calculations. The emit- 
ting surface may be thought of as a number of tiny tips (Fig. 
9) with height of h 2 and sharpness of r 2 . The minimal sepa- 
ration between tips s was assumed to be equal to h 2 because, 
if s<h 2 , a screening effect would take place and lead to a 
decrease in the height of the tips. The maximum number of 
the tips on the surface N is thus approximately equal to 




Fig. 9. Geometry of the emitter, assumed by calculations of emission from 
surface protrusions by using the TSFE approach. 



[(r x /h 2 ) + \] from simple geometrical considerations for 
square r x on the side. Since no sharp surface protrusions 
could be seen at maximum magnification of the SEM, the h 2 
should be less than 20 nm (this is the best resolution in most 
cases of morphological investigations by SEM). 8 The two- 
step field enhancement (TSFE) approach was used in this 
case. The electric field on the blunt tip is equal to 

E x =(h x lr x )E 0 = (h x lr x ){Vld). (3) 

Field enhancement on the very ends of protrusions can be 
written as 

E 2 = {h 2 /r 2 )E x = (h x h 2 /r x r 2 )(V/d). (4) 

In Fig. 8(a) /- V plots calculated by these two approaches 
compared to experimental ones are shown. As can be seen, 
both models could describe experimental behavior. In the 
Table II emission parameters, compared to the calculated 
curves are presented. At present it is difficult to make a de- 
cision as to which of two approaches corresponds to experi- 
mental behavior. For the first approach it is unclear how 
diamond could have a number of electrons in the conduction 
band or near the conduction band sufficient to produce large 
emission currents. For the second approach an explanation is 
needed as to why stable ultrasharp protrusions are formed 
reproducibly on the diamond surface. 

IV. CONCLUSIONS 

Diamond-coated silicon field emitters demonstrate high 
emission currents at moderate voltages by very blunt tips. 
This behavior can be explained as due either to lowering of 
the effective work function or due to a large field enhance- 
ment at the tiny ultrasharp protrusions on the emitting sur- 
face. 
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Engineers are racing to avert what could become a plague of short 
crcuits in electrical and electronic fences. By Ivan ArZato 



in the cold Vacuum, of 
space, on a gleaming metal sur- ; 
face inside the Galaxy 4 com- 
munications satellite, tiny whisk- 
ers of tin grew in perfect stealth 
—until May 19, 1998, that is. 
That's when at least one of 
those whiskers bridged a pair of metal contacts in the 
satellite's control processor. The short circuit killed 
the satellite. Some 40 million pagers stopped work- 
ing all over the country. Millions of dbllars" worth 
of ATM and credit card transactions were inter- 
rupted. The $250 million satellite became, in the 
words of NASA engineer Henning Leidecker, "a ; 
doorstop in space." 

The loss of Galaxy 4 was just one of the more vis- 
ible consequences of a little-understood problem 
with catastrophic potential for electronic and elec- 
trical systems: metal that grows whiskers. An F-15's 
radar system, pacemakers, fuse switches in air-to-air 
missiles, electronic relays in a nuclear power plant 
and global positioning system receivers— not to men- 
tion many other satellites— all have fallen victim to 
the problem. One group of University of Maryland "Z 
theorists has estimated that tin whiskers have caused- " 
losses of billions of dollars to date. 

And the problem suddenly could be poised to get 
worse. The march of miniaturization means that ever 
smaller metal whiskers can short out the ever smaller 
distances between leads, solder bumps, and other 
jam-packed conductive surfaces in electronic sys- 
tems Furthermore, the European Union has or- 
dered thatbythe middle of next year, electronic and 
electrical products sold in its vast market must be 
free of lead— and it is lead .that best checks the- 
growth of whiskers, ^are more vulnerable now 
than we used to be," says Ron Gedney, director- ; 
turned-consultant for the National Electronic Man- 
ufacturing Initiative (NEMI) in Herndon, Va 

The phrase "tin whiskers" is pretty darned de- 
scriptive. Under just the right lighting, if you look at 
a metal surface with whiskers, it sparkles. No one can 
tell you the specific conditions under which whiskers 
may or may not grow, but it usually has something to 
do with surface tension. Justas. you might move to 
the less crowded outside wall of a room packed with 
people, metal atoms move around and reform into 

^^l^^^^^^^press- 
mg of a bracket or screw. 

- Until the new millennium the tin-whisker prob- 
lem actually seemed like a plague of the past The 



mm 1 



ip. 




metalluigical bug first became apparent in the late 
-1940s to telecommunications engineers who were 
investigating why relays in telephone switching sys- 
tems were failing. The research community came up 
with a whisker-quashing solution: Add 2% to 3% 
lead to the tin plating used in electronic assem- 
blies, particularly on wires and leads to make them 
solderable. Lead-tin aUoys became standards for the 
mdustry, and their use relegated tin whiskers, at least 
tor most of us, to nuisance status. 

, Asoftenhappens,thesolutionbecameaproblem- 
.Lead became amaterial/Kw^naB and for good rea- 
; son. It does a nasty number on neural circuitry es- 
pecially if you're a kid. To reduce human exposure 
to lead, governments have regulated it out of paint 
■ gasoline, plumbing, and other sources. The relatively 
- small amount of lead melec^mc and electrical sy/ 
. tems-^about 0.2 gram in an Intel processor and 
about two to three grams in a motherboard-hadn't 
attracted much attention, but that changed as the vol- 
ume of electronic waste showing up in landfills be- ' 
gan balloonmg. 

Although some companies had been anticipat 



Metal whiskers grow 
In many shapes. The 
most common is the 
"needle" (top toft), i 
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FIRST 



this Just In 



ing an era of lead-free electronics for more 
than a decade and already are shipping 
lead-free products, the get-the-lead-out 
clock started ticking in earnest for the en- 
tire electronics industry in 2002 when the 
European Union enacted the Restriction of 
Hazardous Substances (RoHS) and the 
Waste Electrical and Electronic Equipment 
(WEEE) directives. Among other things, 
me directives mandate that by July 1, 2006, 
only lead-free electronic and electrical 
products will be legally salable 
in member countries. Unfortu- 
nately, no one knows for sure 
whether there is a no-lead sub- 
stitute that will ;be immune 
to whisker growth. And that 
means tin whiskers and the 
failures they can foment could ; 
become as familiar in elec- 

; . tronic products as flat tires are 

'.'in cars. '■'};•■'' 

To those who make ele:c^ 
trcmic products, the no-lead re^; 
striction is akin to phasing salt 
out of food. Lead is a trace in- 

: gradient in these products, but 
its whisker-stbp>ping and other traits have 
a technical taste that's hard>tp. do without. 

^ Even though some comply 
ing to have developed rehabM substitutes 
to lead-tui ialld^s; a fundamental metallur- 
gical ignorance about why whiskers form 



Getting the 
lead out of 
electronics 
may aggravate 
the whisker 
problem. 



even a decade ago, well before the heavy 
metal made it to the top of regulators' 
to-do lists. By the late 1990s the likelihood 
that lead would be excised from the diet 
of electronics makers and suppliers be- 
came more apparent, and in 2002 the Eu- 
ropean Union made it official. Japan, 
China, and several American states, in- 
cluding California and New Jersey, also 
are in various stages of outlawing lead. 
Denial is not an option. 

Almost as out of sight to most 
people as the tin whiskers them- 
selves is a worldwide network of 
detail-minded engineers, met- 
allurgists, technology managers, 
chip manufacturers, govern- 
ment officials, and others who 
are racing against the clock to 
make sure that the get-the- 
lead-out movement ddes not 
awaken the tin-whisker driigcm. 
They have until about right 
now to come up with a solution 
if they want to make sure that 
the products^that ship in the 
coming months and years do 
not collectively harbor a Y2K of sorts-^the 
wide distributicm of everything horn mi- 
crowave ovensUo missiles that are^more 
prone.tq whisker-induced failure than most 
eiectronicproducts have been for the past 
hatf-century:<v "_„, irr - - /- "■■ 




on metd^mcluding cadmjum; zinc, and >, "You don't wantHb^go. out and buy a 
.silver in addition to tin---means that there 'brand-new HDTV, only for it to fail before 
is ho rock-solid basis ; for expecting success you've finished your payments," says Ged- 
through any particular fix. That's why ney. JResearchers have come up with a va- 
military:systems are, for the moment, ex- j riety of potential solutions. Semioonductor 
_ empt from the lead-free requirements. But ^ maker, Ag^ 

' "becau^ military procurement has relied nduhced in September that it will under^ 
more and more on off-the-shelf porTsum^r coat the leads on the components it uses 
products; ^l^iffmilitary systems'in coming with|nkkel before it puts pure tin on top. 
years could end up with more vulnerable -- yEhgjm a tin-siiver- 

f comF^nents. ^^ copper alloy that appears to limit whiskers 



"Wl^np^ha^^^ to amo^mamag^ble fact of metallic life, 

and you have 50 years of experience with 
it, yot hate having to change it," says 
NEl^'s^^^What's twbad, he' and;. 



eFpromising those iead-free so- 
lutions might look in laboratory and beta 
tests, however^the reallssessment of their 
long-termtability to keep whiskering at bay 
will be taking place oveyhe riext few years 



others-in the mdiistry a%ue, c is that going 
after lead in electronics is not a good way 
^bffgetting a lot of environmental bang for^in^the wild asj^d^e electronic systems 
the buck^especially when it includes the* are madefahia sold throughout the world, 
riskof brmging down more satellites. Even> It^s#Sone deal for any electronics indus- 
sorSays Gedne^no^fi^^ that wants to remam in the game, 

against progress on the env^nmentaf ^Init all are moving forward with the nag- 
^ront^ei^ ging anxiety that they might be setting tril- 



a^soiutiorijto thisf \ s- lions of individual stages for the quiet, 

/ Some / comjpanies, such as Texas Instru- , stealt^grdwth of metal whiskers that can 
mentsf developed lead-free technblpgies ^^o no good. Q 
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Physics 114 Spring 2000 Week 6.2 
Muons in Matter 

You just determined how the binding energy of an electron to a unit positive charge is changed 
by the dielectric constant of the material it's in. That problem is crucial for understanding how 
semiconductor devices work, Here's another problem that's important for understanding how 
certain types of analysis of materials work In these techniques, particles called muons are 
injected into materials, and the types of sites that they bind to are studied. 

A muon is like an electron but with 200 times as much mass. 

(Hint- you may need to refer to the algebraic solution for the radius, etc of a hydrogen-like 
atom in terms of the mass, charge, ...) 

a) What is the radius of the ground state of a muonic hydrogen atom, in which the muon 
replaces the electron? Express your answer in terms of Bohr radii. (Don't worry now about the 
fact that the muon mass isn't negligible compared to the proton mass.) 

As we found in part a ) of problem 1 , 

h 2 

a cc — — 
he m 

Between the electron and the muon, the only quantity in a that differs is the mass, m. Thus, 

a electron *W 200 

a muon = Electron/ 200 = a 0 ' 200 

Note that this uses the correct value for a electron = aQ, not a dectron = 2a Q , the value found in problem 1 
using the approximation to the energy of the hydrogen atom ground state. 

Also note that the advice given above not to worry about the fact that the muon mass isn't negligible 
compared to the proton mass refers to the fact that the m of which we speak is in fact not the mass of 
any single particle but rather the reduced mass of the atom. Using the reduced mass accounts for the 
fact that, considering the issue in classical terms, the electron does not simply orbit a stationary 
nucleus, but rather the electron and nucleus orbit their common center of mass. For atoms containing 
electrons, the reduced mass is very nearly the same as the electron mass, so we ignore this issue 

b) What is the binding energy of a muon to a proton, in eV? 

As we found in part b) of problem L 



http://wug. physics.uiuc.edu/courses/phys 1 1 4/summer02/Discussions/html/wk6/ sol6_2. htm 
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— — 
h 

Between the electron and the muon, the only quantity in E that differs is the mass, m. Thus, 

Eekctren m electron ^ 

The binding energy of a particle is defined as the minimum energy that must be supplied to free the 
particle. All bound states of atoms have negative energy, so the muon is free as long as its energy is 
greater than or equal to zero. Thus, the binding energy of the muon, assuming that the muonic 
hydrogen is in its ground state, is the absolute value of the ground state energy. 

Ending, muon = l E muo„l = 200 l E electn„J = 200 H 3 6 cV| =2720 eV 

This uses the correct value for the ground state energy of the hydrogen atom, -13.6 eV, not the value 
found in problem 1 using the approximate expressions. 

Note that when I write of the energy of the electron or muon, I really mean the energy of the hydrogen 
or muonic hydrogen, because as was mentioned above, not only the electron or the muon but also the 
proton has momentum. 
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EFFECT OF SCREENING ON THERMONUCLEAR FUSION IN STELLAR AND LABORATORY PLASMAS 

L. Wblets 1 , B. G. Giraud 2 , M. J. Watrous 2 , and J. J. Rehr 1 
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ABSTRACT 

The fusion-enhancement factor due to screening in the solar plasma is calculated. We use the finite- 
temperature Green's function method, and a self-consistent mean field approximation. We reduce this to 
one-center problems, because in the collision of two fusing ions, the turning point where tunneling may 
occur lies far inside the screening radius. The numerical results given by this method indicate that 
screening may be slightly weaker than that obtained in the most recent previous calculations. 
Subject headings: elementary particles — nuclear reactions, nucleosynthesis, abundances — plasmas 



1. INTRODUCTION 

There remains continued interest in the solar neutrino 
problem, which has not been resolved in terms of the stan- 
dard solar model and standard particle physics. Neutrinos 
are generated in thermonuclear reactions, and thermonu- 
clear reaction rates are determined by solar physics 
(temperature and density), nuclear physics, and atomic 
shielding. The solar neutrino problem is just one of several 
examples where there is current interest in thermonuclear 
reaction rates; others include general stellar structure and 
laboratory plasmas. We address here the role of atomic 
shielding, using a finite-temperature Green's function 
method. Recent calculations of screening effects in the Sun 
have been reported by Gruzinov & Bahcall (1998), who also 
review previous solar screening calculations. Results of 
those calculations, with which we compare our present 
work, exhibit a scatter of 10%-30% or more in the effect of 
screening on fusion rates. As discussed in Gruzinov & 
Bahcall (1998), all of these previous calculations either are 
based on Salpeter's weak screening formula, which is a clas- 
sical approximation, or else make unphysical approx- 
imations in order to go beyond the linear regime. Gruzinov 
& Bahcall attempted to rectify these limitations by includ- 
ing quantum effects, treating the kinetic energy in pertur- 
bation theory. Our model differs from theirs in that we do 
not rely on perturbation theory. 

The physical model we employ is based on self-consistent 
Hartree plus finite-temperature local density approx- 
imation for exchange and correlation, and it assumes 
thermal equilibrium and adiabaticity in the internuclear 
coordinates. We solve this model making no further 
approximations or expansions, except for numerics that are 
carefully monitored. AH previous calculations involve, e.g., 
expansions in the temperature, density, or quantum correc- 
tions. Carraro et al. (1988) also consider dynamic effects. 
Thus we believe that our approach yields results that are 
the most accurate to date. These results yield fusion- 
enhancement factors that are typically a little smaller than 
those recently reported (Gruzinov & Bahcall 1998). 

The rate of thermonuclear fusion in plasmas is governed 
by barrier penetration. The barrier itself is dominated by 
the Coulomb repulsion of the fusing nuclei Because the 
barrier potential appears in the exponent of the Gamow 
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formula, the result is very sensitive to the effects of screening 
by electrons and positive ions in the plasma. Screening 
lowers the barrier and thus enhances the fusion rate; the 
greater the nuclear charges, the more important it becomes, 
and thus it plays an important role in the solar neutrino 
spectrum. 

As shown below, the classical turning point radius that 
enters in the WKB integral for barrier penetration is very 
small compared with the characteristic screening lengths of 
interest. Inside this radius, the barrier potential is just the 
nucleus-nucleus Coulomb repulsion minus a constant due 
to screening. The constant can be interpreted simply as the 
difference in free energies of the system between united and 
separated nuclei. Since both of these are spherically sym- 
metric, one needs only consider one-center problems. 

We present here numerical calculations relevant to the 
solar core. The screening due to electrons in the plasma is 
calculated quantum mechanically by a novel Green's func- 
tion method described in a recent paper (Watrous et al. 
1999). The screening due to the ions is treated classically but 
self-consistently with the electrons. The resultant enhance- 
ment factors for several nuclear reactions are presented and 
compared with earlier results by other authors. 

2. FUSION RATE 

We work in atomic units, e = h = m e = k B = 1 ; the unit 
of atomic temperature is 3.159 x 10 s K. As shown by Salpe- 
ter (1954; here we follow Clayton 1968, p. 359), the fusion 
reaction rate between species 1 and 2 is given by 

r U2 = N t N 2 ^WMmEWE , (1) 

where N t and N 2 are the number densities of the colliding 
nuclides, #E) is the Maxwellian probability that the center- 
of mass energy at large separation is £, and the cross section 
a is written as a product of a penetrability factor P and a 
nuclear factor <r nuc , 

a(E) = P(E)a BUC (E). (2) 

In the WKB approximation, the s-wave penetration is given 
by 

x exp [-ZJ% J^[Z l Z 2 /r + UJr) - KJ^drJ, 

(3) 
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with E B the height of the barrier, \i the reduced mass, and 
tfsc(r)<0 the screening potential; R is the "touching" 
radius (the top of the barrierX and R 0 is the classical turning 
radius. The integrand in equation (1) peaks at some charac- 
teristic energy F 0 . For the solar interior, the corresponding 
classical turning radius R 0 « Z i Z 2 /E 0J is ~ 0.005 AU, 
compared with a screening radius of ~0.5 AU. For this 
reason, the screening potential can be taken as a constant 
U 0 = U sc (r = 0) < 0 within the turning radius, acting as an 
effective shift in the energy : 

r ta = N t N 2 £m*mE-U 0 )*^E-U 0 )dE (4) 
r u2 oc £V/ 2 <T E/r P(£- U 0 )c nuQ (E-U Q )dE (5) 
= P° (E + UJ l *e-* + WP(BI*^tBWB • («) 

J-Uo 

Then, neglecting U 0 in the limit of integration, 

r lt2 a e- v « jT N t N 2 J" " \KEHF)a(E)dE . (7) 

Screening thus enhances the reaction rate by the factor 

/=e- Uo/r , (8) 

where U 0 is given by equation (10) below. Since U 0 < 0, 
/ > 1, corresponding to enhancement. 

3. FREE ENERGY AND THE ONE-CENTER PROBLEM 

The screening potential is just the change, brought about 
by the approaching ions, in the Helmholtz free energy 

F = Internal " TS , (9) 

as a function of nuclear separation. Because of the smaDness 
of the turning radius, it is sufficient to consider the one- 
center problems and identify 

U 0 - F{Z V + 2 2 ) - F(Z t ) - F(Z 2 ) . (10) 

What we mean here by one-center problems are the calcu- 
lations of F when either charges Z ls Z 2 , or both lie at the 
center of the plasma. In § 4 we calculate F(Z) explicitly for 
an independent particle model approximation, akin to the 
temperature dependent Hartree-Fock method. However, 
the free energy has a simple and completely general (for the 
problem at hand) interpretation, which does not depend on 
the mean field approximation. It is, by definition, the work 
required at a given temperature to introduce the nuclear 
charge into the plasma. The nucleus interacts with the 
plasma only electrostatically, so dF = #(Z', r = 0)dZ', and 
thus the work required to assemble the charge Z at, say, 
r = 0is 



F(Z) 



-f 



4>(Z',r = 0)dZ, 



(11) 



where <f>(Z\ r) is the electrostatic potential generated in the 
plasma due to a nuclear charge Z', 

4>{Z\ r = 0) - An J["**[z - P&> d] • < 12 ) 

Here pjr, T) and pfa T) denote the finite-temperature 
number density (recall e = 1) of the electrons and various 
background ions charge densities, respectively. 

The one-center screening problem is solved by the 
methods described by Watrous et aL (1999) for the electrons 



with a Kohn-Sham formalism, using the finite-temperature 
local density approximation (Mermin 1965) for exchange 
and correlatioa Screening due to ions is included as well as 
that due to electrons. The ions are treated classically 
according to the Debye-Huckel method. Thus the electric 
potential generated by electrons and ions becomes 



am Z _, f , jy l/P/^r)-^, T) 



(13) 



The second term is the function denoted by <£(Z, r) above. 
Each ion density is given by 

pAr,T) = pAoo)e- z <™ T . (14) 

Since tyr) is positive, ions are pushed outward and electrons 
are drawn toward the nucleus. For notational simplicity in 
the following, the temperature dependence of the densities 
will be understood and not specified explicitly. 

4. ALTERNATIVE MEAN FIELD DERIVATION OF THE 
FREE ENERGY 

We now demonstrate the equivalence of the F defined by 
equations (9) and (11) in the finite-temperature Hartree- 
Fock approximation, which is similar to the finite- 
temperature local density approximation actually used in 
our numerical calculations. In a transparent, short nota- 
tion, we describe the system by a second-quantized Hamil- 
tonian containing the kinetic energy operators t t of the 
particles present in the plasma, their two-body interactions 
v ij9 and the additional one-body operators = — Z/rj or 
-hZjZ/r,- representing the contributions of the additional 
nuclear charge Z at the center. The temperature T enters 
the formalism via the usual Boltzmann factor exp 
where P = 1/T. Actually, because the average electronic 
density and background positive charge density are fixed 
parameters of the problem, Jf must be replaced by the 
usual constrained Hamiltonian. This will read, in a short 
notation, — \j!JT. Here p! is a several-component chemi- 
cal potential because of the several components in the 
plasma (ie., electrons and various background ions). 
Accordingly the particle number operator Jf must be 
understood as a several-component operator. For simpli- 
city, however, the following equations deal with the elec- 
trons only. This is because the positive ions in the 
background may consist of bosons as well as fermions, and 
we want to spare the reader the cumbersome symmetriza- 
tion or antisymmetrization formulae for such backgrounds. 
In any case, the present paper describes the background 
density classically, and this density is not high enough to 
demand exchange terms. 

The independent particle Ansatz approximates the eigen- 
states of — fi'^V by antisymmetrized products of single 
particle states (orbitals) | i>» with an approximate spectrum 
made of sums of independent particle levels t\i. The infinite 
set of orbitals | i> and eigenvalues % defines a one-body 
operator Jf 0 — h'jV, presumed optimal from the point of 
view of a variational principle: Le., stationarity of the grand 
potential ot free energy F. 

In second quantization, the many-body density matrix 
used as a trial density operator (noted in the following by 
quantities carrying a subscript t) is 

exp l-P(J? 0 - liJT)] exp frcfc) 

D t - z z ' * 
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TABLE 1 

Results of Solar Screening Calculations for 
Various Values of Z 

Z <ft(Z,r = 0) pJjZ,r = 0) F(Z) 

1 -2.007 2.791 -1.008 

2 -3.944 6.829 - 4.009 

3 ... ... -8.962 

4 -7.884 20.76 -15.82 

5 ... ... -24.54 

6...... -1Z00 49.77 -35.14 

7 ... ... -47.65 

8 -16.51 108.6 -6219 

Note. — <f> and p e are calculated values; the free 
energy F is based on the polynomial fit 



where c\ and c t are the usual fermionic creation and anni- 
hilihation operators for orbital 1 1>- For electrons the nor- 
malization denominator Z t is 



Z, = n[l + e*P(-M], 



so the trial density becomes 



n exp(-jgfrcfo) 
"'"M l + exp*-/?*) 



(16) 



(17) 



The true density matrix D = exp - n'-^)]/ 

Tr exp [-/?(Jf - tfJf)\ minimizes the true grand poten- 
tial, n = Tr D(Jf - v'JT + J9~ MnD). However, with an 
independent particle approximation, one must be content 
with minimizing 

F = Tr Dpf - n'JT + /T 'InD,) . (18) 

It is easy to show that 

Tr D t c}c iC }cj f p if i ¥>j , (19) 

where the Fermi occupation numbers are^ = Tr Djc/c, = 
1/[1 + exp (prii)]. Furthermore 

Tr D t {-faicjci - ln[l + exp (-/&&)]} 

^ta/iH- (l-^lna-ZJ. (20) 

Finally, because = + t^*) + ^ is the sum of a one- 
body operator 2T + if = t + I,- w 4 and a two-body 
operator Y = £,•>,• we obtain 

f = Z /;•<*!(* + w - joio + i z y; /xoi^iy> 
+ r *i UiWi + a - ami -jm • (2D 

i 

Here the tilde t; means that the matrix element <(/ | v | y> is 
antisymmetrized. 
The functional derivative of F with respect to <i | is then 

=/i (t + w - u' + £ /^j|i>|y>) , (22) 

where one recognizes the action of a mean field potential, 
including both direct and exchange terms 



Stationarity of F with respect to |i> then gives the "finite- 
temperature Hartree-Fock" equations, 



(t + w - fx' + <&0|i> = VilO. Vi > 



(24) 



where one recognizes that r\i is the Lagrange multiplier for 
the normalization of the orbital 

In the same way, the derivative of F with respect to ij h or 
as well f h yields 



^ = ^ i + r l inc/ i /i -/J, 

ok 



(25) 



which vanishes [see eq. (20)]. It can thus be concluded that 

F = -/T 1 InZ, - 1xD t ir , (26) 

(which differs from - P" l h\Z t ) is stationary with respect 
to variations of D v Therefore, since if is proportional to 
the additional charge Z at the center, the derivative of F 
with respect to Z must be given by 



dF ^ ^iT 
3Z = TTD <Y 



(27) 



Notice also that iT is a local potential Hence, in the coor- 
dinate representation, only the diagonal part of D t (the clas- 
sical density) is needed to calculate dF/dZ, 



|| = 4^Jr'^[p p (/)^p e (rO], 



(28) 



W\i> = Zfj<j\Z\ij> 
j 



(23) 



which is the differential form of equations (11-12). Here 
pjj*) and pjj-') are diagonal matrix elements <r' | D t | r'> in 
the electron and positive background sectors, respectively. 
It will be noticed that equation (28) gives the background 
density as well as the electronic density, while the preceding 
equations, equations (15H27), accounted for the electrons 
only. In view of the simplicity of the result due to the elec- 
trons, this reinstatement of the background contribution is 
trivial. Notice also that for the one-center problem, D t is 
rotationally invariant, hence no vector label r' is needed. 

The result, equation (28) is simply a reformulation of 
equations (11M12). The method of Matsubara poles used 
by Watrous et al (1999) via the finite-temperature, one- 
body Green's function, is perfectly suited to this local calcu- 
lation of the density and avoids an explicit solution of 
equation (24). Indeed, as discussed in detail by Watrous 
et al. (1999), a local density such as p e (r') can be directly 
derived from diagonal matrix elements 
<r' | (i; - t - w - p' - <%y 1 1 r'> of the one-body Green's 
function. Such matrix elements are integrated on a suitable 
contour in the fj-complex plane. Once dF/dZ is known for 
all values of Z smaller than Z t + Z 2 , trivial integrals 
provide the screening U Q according to equation (10). 

5. NUMERICAL RESULTS FOR THE SOLAR CORE 

Shielding calculations have been performed for Z - 1, 2, 
4, 6, 8 at a density and temperature relevant to energy and 
neutrino production in the Sua From these calculations the 
fusion-enhancement factors are calculated and compared 
with the results of other researchers. Ricci et al. (1995) have 
shown that enhancement factors for the relevant fusion 
reactions are insensitive to the location within the solar 
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TABLE 2 



Fusion-Enhancement Factors for the Solar Interior (This Work), 
Compared with Other Screening Calculations 



Reaction 


This Work 


WES 


Mit 


GLXxC 


CSK 


GB 


P + P 


1.043 


1.049 


1.045 


1.049 


1.038 


1.053 


He + He 


1.181 


L213 


1.176 


1.115 


1.158 


1.224 


Be+p 


1.183 


1.213 


1.171 


1.112 


1.169 


1.166 


N+p 


1356 


1.402 


1.293 


1.191 


1.324 


1.393 



Note. — WES is Salpeter's weak screening approximation (Salpeter 1954); Mit is 
Mitler's Thomas- Fermi-like model (Mitler 1977); GDGC is from Graboske et aL 1973; 
CSK is the dynamic screening model of Carraro et aL 1988; GB is from Gruzinov & 
Bahcall 1998. 



core. We use here the set of parameters similar to those 
employed by Bahcall & Pinsonneault (1995), corresponding 
to R/R 0 = 0.06. In atomic units, they are T = 47, n e = 7.63, 
X = 0.432, Y = 0.568. The results are displayed in Table 1 
for 0(Z, r = 0) and Ap^Z, r = 0) = p«(Z, r = 0) - n c . These 
quantities for other Z-values can be obtained by inter- 
polation. 0(Z, r = 0) was fitted to a fourth-order poly- 
nomial, 

#Z,r«0)«-£c ll Z". (29) 
«=i 

The integral of this quantity then yields the free energy, 

f(Z)*-£*Z» +1 A»+l). (30) 

with c t = 2.0275, c 2 = -0.03661, c 3 = 0.00594, c 4 = 
- 0.000103. Note that 0(Z, r = 0) is nearly linear in Z. 

Table 2 gives the fusion-enhancement factors for several 
reactions of interest in the solar neutrino problem. The 
deviation from unity is a result of shielding. Note that the 
various calculations exhibit a scatter of 10%-30% and 
more in the deviation from unity. Although our factors are 
rather close to the recent calculations of Gruzinov & 
Bahcall (1998), they are usually somewhat lower. 



6. SUMMARY AND CONCLUSIONS 

The finite-temperature Green's function method of 
Watrous et aL (1999) has been applied to the problem of 
screening of the nucleus-nucleus interaction in the solar 
plasma. The method is based on a self-consistent, finite- 
temperature Kohn-Sham formalism with the local density 
approximation for exchange and correlation. Atomic bound 
states are included on an equal footing with continuum 
states. Fusion-enhancement factors (over pure Coulomb) 
are calculated for various relevant nuclear reactions at 
mean conditions (temperature and density) of the solar core. 
Comparisons with several other calculations are presented. 

The method appears to have no restriction with respect 
to temperature or density for stars in "noncatastrophic," 
thermal equilibrium states. When implementing calcu- 
lations for other systems, one may need to extend the range 
of parameters for the local exchange-correlation term, 
which is quite small for our solar system. 

We wish to thank Professor G. Krein for valuable dis- 
cussions on all aspects of the paper. 
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